Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 16 November 2009 (MN WI^ style file v2.2) 



Reconstructing the Arches I: Constraining the Initial 
Conditions 



o 
o 

> 

o 



a 

o 



> 

oo 
in 
o 
m 



o 



S. Harfst^* and S. Portegies Zwart^ and A. Stolte^ 

^ Sterrewacht Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands 

^1. Physikalisches Institut, University Cologne, Ziilpicher Str. 77, 50937 Koln, Germany 



16 November 2009 



ABSTRACT 

We have performed a series of A^-body simulations to model the Arches cluster. Our 
aim is to find the best fitting model for the Arches cluster by comparing our simulations 
with observational data and to constrain the parameters for the initial conditions of 
the cluster. By neglecting the Galactic potential and stellar evolution, we are able to 
efficiently search through a large parameter space to determine e.g. the IMF, size, and 
mass of the cluster. We find, that the cluster's observed present-day mass function 
can be well explained with an initial Salpeter IMF. The lower mass-limit of the IMF 
cannot be well constrained from our models. In our best models, the total mass and 
the virial radius of the cluster are initially (5.1 ± 0.8) • 10^ Mq and 0.76 ± 0.12 pc, 
respectively. The concentration parameter of the initial King model is wq = 3 — 5. 

Key words: stars: formation - stellar dynamics - methods: A^-body simulations 



1 INTRODUCTION 

The Arches cluster is one of only a few young and mas- 
sive starburst clusters in the Milky Way. Its location at 
a projected distance of less than 30 pc from the Galactic 
center and an age of only ~ 2.5 Myr ijPiger et al.l l2002l : 



iNaiarro et al.|[2004l 'l make this cluster a unique object for 
studying star formation and dynamical processes in the cen- 
ter of galaxies. 

The observed present-day mass of t he Arches cluster 
' has been estimated w ith ~ 1 - 2 ■ 10'* M© (jPieer et al.ll 19991 : 

lEspinoza et al.ll2009l ). With this mass a cluster will not sur- 
vive long in the Galactic center environme nt and evaporate 
on a time scale maybe as fas t as ~ 10 Myr ()Kim et al 1 l 19991 : 

IPortegies Zwart et al.|[2002l ). The initial mass of the clus- 
ter has been determined from A'^-body simulations, how- 
ever, different results h ave been obtained by different au- 
thors: iKim et all (|2000| ) found that their best model for 
the Arches cluster had a total mass of about 2 • 10* Mq; 

IPortegies Zwart et al.l l|2002h . on the other hand, came to 
the conclusion that the cluster was initially more massive 
than ~4 - 1O'*M0. 

The initial mass function (IMF), a key aspect of star 
formation, see ms to be uniform throughout the universe 
l|Kroupa|[2003 ). This universal I MF c an be described by 
the power-law found by Salpcte for stars in the so- 

lar neighborhood and is valid from 0.5 Mq to the highest 
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masses. Below 0.5 Mq, the IMF is significantly flattened (e.g. 
iKroupalliool l. 

Determining the IMF of young clusters from observa- 
tions is not a straight-forward process. Uncertainties can 
arise from the measured luminosities, the estimated age of 
the cluster, the completeness of the observed sample, and 
the stellar evolution models. In addition, the non-linear dy- 
namical evolution of the cluster has to be taken into account 
as shown in Fig. [1] as the star cluster evolves, more massive 
stars (star symbols) will move towards the cluster center and 
low-mass stars (points) in the opposite direction (indicated 
by the arrows in the left image). If the detection of clus- 
ter members is radially limited (dashed circle) , it will result 
in an observed mass function (MP) in the mass-segregated 
cluster (right image) that is different from the IMF. This 
effect is visualized in Fig. [1] by the ratio of low- to high- 
mass stars inside the dashes circles before and after mass 
segregation. 

In case of the Arches cluster, observations have revealed 
that the slope of the observed mass function (MP) is sig- 
nificantly fiattened with P ~ —0.9 ± 0. 15 with respect to 
the standard Salpeter IMF (P = -1.35) (|Stolte et al.ll2005l . 
l2002l : |Piger et al.^ll999^ ■ and therefore the Arches cluster has 
been regarded as an argument against the universality of the 
IMP. More recently, however, Espi noza ot a l. (2009) derived 
a slope of P = —1.1 ± 0.2 in i? < 0.4 pc and concluded that 
a standard Salpeter IMP cannot be ruled out for Arches. In 
addition to the radial variation in Av, these authors also 
accounted for differential extinction variations, which can 
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Figure 1. Schematical view on the cluster mass function and 
evolution. Low- and high-mass stars are shown by points and 
star symbols, respectively. The dashed circle indicates an radial 
observational limit. In the left image, the cluster is at t = 
with arrows denoting the effects of dynamical evolution. The right 
images shows the mass-segregated cluster. 



severely affect the incompleteness and may have biased the 
earlier results. Large uncertainties in the slope still remain, 
revealing the necessity to compare the observed cluster MF 
with simulations. 

In addition to the flattened slope, there has been some 
debate whether the IMF of Arches is truncated at the low- 
mass end as the result of the extreme conditions at the 
Galactic center where the clus ter has form e d. Po ssible evi- 
dence for this was reported bv lStohe et aH (|2005| ). who de- 
termined a low- and intermediate-mass depleted MF in the 
cluster core (R < 0.2 pc) with a turn-ove r at 6 — 7 Mp) . Thi s 
truncation in the MF was not seen by iKim et al.l 1 200d ). 
They only found a local bump in the MF at ~ 6M0. Even 
if the MF is truncated at the low-mass end, it remains un- 
clear whether this would be the result of a truncated IMF or 
a dynamical effect such as tidal stripping of low-mass stars. 

Several studies have been done to determine the IMF of 
the Arches cluster using numerical simulations, again com- 
ing to different conclusions: the model favored bv lKim et al.l 
l|2000t ) started with fla t IMF with a slope of F = -0 75 close 
to the observed one. iPort egies Zw art et al.l (|2002| ) found 
that the observed IMF is the result of observational selection 
effects and dynamical evolution of the cluster. They argue 
that the observed flat MF is therefore co nsistent with a nor- 
mal Salpeter IMF. And lKim et all l|2006l ) suggested that the 
initial slope of the MF was —1 to —1.1, only slightly shal- 
lower than Salpeter. 

The Arches cluster also exhibits clear signs of mass 
segregation. The slope of the observed MF for stars 
in different annuli changes with distance from the clus- 
ter center. Towards the center, the slope becomes shal- 
lower and further out the slope is closer to Salpete r 
llKim et al.ll2006l: IStolte et al.|[2005l : lEspinoza et al.ll2009l ). 
IPortegies Zwart et al. I l|2007h have found the same charac- 
teristics in numerical A'^-body models and concluded that 
the Arches cluster is about half-way to core collapse. 

Despite the wealth of detailed observational data, large 
uncertainties regarding some of the properties of the Arches 
cluster remain, most importantly in the slope of the observed 
MF. Numerical simulations have been used to better under- 
stand the observations, but no fiducial model has emerged 
from these studies so far. In this paper, we want to recon- 
struct the initial properties of the Arches cluster by com- 
paring the results of A'^-body simulations with observations. 
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Figure 2. The completeness fraction in two radial bins as func- 
tion of initial stellar mass. The data points (blue c rosses and red 
pluses) are from the analysis of lStolte et al.l ll2005h . The two full 
lines show a fit to the data. 



In our A''-body simulations, we model the Arches cluster 
on a star-by-star basis and, for the first time, systematically 
explore the parameter space to find the best set of initial 
conditions describing the Arches cluster. We vary parame- 
ters determining the initial mass, size, and concentration of 
the cluster. In addition, we also test which IMF can best 
explain the MF observed today. The total number of free 
parameters in our models is five, and a large number of 
simulations is required to search the full parameter space. 
Therefore, we decided to neglect the Galactic potential and 
stellar evolution for now. Once the initial conditions of the 
cluster have been constrained, we will use these results for 
simulation that include these processes in a follow-up paper. 

The paper is outlined as follows: in Sec. [51 we describe 
the observational data we use. Then we explain the param- 
eters for the cluster model in Sec. |3] and the simulations and 
their results in 3] We conclude and summarize in Sec.[S]and 
Sec.El 



2 OBSERVATIONAL DATA 



We use data from observations bv lStolte et al] l|2005l . NACO 
data hereafter). The NACO data has been taken using the 
ESO VLT AO system NAOS and the CONICA near-infrared 
camera in two wave bands, H and Ks- The field of view is 
~ 25" squared or 1 pc^ (we assume a distance of 8 kpc to 
the cluster) with a resolution of ~ 0."08. 

The total data set consists of ~ 2200 stars belonging to 
the cluster and the field. In order to fi nd the cluste r stars, 
we apply the same colour selection as IStolte et al.1 (|2005| ) 
(see their Fig. 2), leaving ~ 1500 stars. We then fit the K- 
band magnitu des (corrected for t he observed radial variation 
in extinction (jStolte et al. 20021)') against a 2.5 Myr Geneva 
main-sequence isochrone (ILeieune fc Schaererll2001 ], assum- 
ing solar metallicity for Arches i Naiarro et al.ll2004 ). From 
this we get both present-day and initial stellar masses for 
each star in the sample. 

Because the incompleteness of the data due to crowding 
effects increases significantly for stars towards the low-mass 
end we reduced the sample further to about 300 stars by 
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Figure 3. The cumulative mass profile for stars with m > 10 Mq . 

selecting only massive stars with m > 10 M0. With our 
chosen lower mass limit, the data is 80% complete or better 
for any star mass and position, as a detailed analysis by 
IStolte et al.l (|2005h has shown. Fig.[2]shows the completeness 
fraction as a function of initial stellar m ass for two di f ferent 
radial bins. We have fitted the results of IStolte et al. I l|2005l . 
blue crosses and red pluses in Fig. [2]) and also extrapolated 
for stellar masses above ~ 15 M0. In the following analysis, 
we use this information to either correct the observational 
data or by randomly removing stars from our models. Note 
that only a few stars are added or removed by this correction 
for m > IOM0. 

The field-of-view of the NACO data is such that only 
within a limited radius all stars of the cluster can be seen. 
This radius is 0.4 pc (see left panel in Fig. 1111 which shows 
an image of all the ~ 1500 cluster stars; the radius of 0.4 pc 
is indicated by a circle centered on the center of density). 
We therefore also limit our sample to the 234 massive stars 
within this radius. The cumulative mass profile for these 
stars is shown in Fig.O where we also show the mass profiles 
for three different King models (iKindl 19661 ). All models have 
the same virial radius rvir = 0.5 pc and are scaled in mass 
to approximately fit the observed profile. The present-day 
profile of the Arches cluster is best described by a King 
model with Wo ~ 7. 



3 THE MODEL OF THE ARCHES CLUSTER 

In order to construct a model for the Arches cluster, we 
compare the results of A'^-body simulations with the obser- 
vational data described above. The simulations start from a 
set of initial conditions with a number of parameters that 
can be varied to find the best fitting model. We have chosen 
the IMF, total mass, concentration, and size of the cluster as 
the free parameters. Other parameters are fixed: we assume 
an age of 2.5 Myr and a distance of 8kpc to the cluster. 

3.1 Initial Mass Function 

It has been discussed whether the IMF of the Arches cluster 
deviates from the norm, with a flattened distribution (for 
massive stars) with respect to a Salpeter IMF. Determining 
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Figure 4. The cumulative observed MF of massive stars (m > 
10 Mq) in the Arches cluster (red dots). The dotted blue line is a 
fitted power -law MF and the g reen dashed line is a Salpeter MF. 
Data is from lStolte et ahl ||2005|) and corrected for incompleteness. 



the IMF of a star cluster is not a straight-forward process. 
The present-day mass function (PDMF) of a cluster is the 
result of stellar evolution and dynamical effects. In addition, 
selection effects then shape the observed MF, which will 
use in the following for the comparison with our models. 
We do not take into account stellar evolution, and we use 
the initial masses of stars to determine the observed mass 
function (using the present-day masses of stars would change 
the MF only for stars with m > 50 Mq). We also correct the 
observed MF for the incompleteness of the data, however, 
the observed MF can differ from the underlying IMF due to 
the selection of stars inside R < 0.4 pc 

We derive the observed MF from the NACO data and 
also find a shallow slope with F — —0.9 ± 0.15 for stars with 
m > 10 Mq in Fig. |4] We correct for the incompleteness of 
the data by weighting stellar masses with the inverse of the 
completeness fraction given in Fig. [2] The slope is derived by 
fitting an integrated power-law MF with two free parameters 
(normalization and slope) to the cumulative observed MF. 
This allows a more straight-forward fitting of the data that 
the commonly used mass binning. The formal uncertainty 
of the fit is small (see Fig. U)) and not tak ing into account 
errors in determining stellar masses. IStoltd (|2003 ) has given 
an estimate of ±0.15 for the total error, which we also use 
here. 

The slope we find here is in agr eement with what has 
been reported bv IStolte et al.l (|2005l ), and significantly de- 
viates from a standard Salpeter MF. As a test, we have 
created random realizations of the Salpeter MF using 1 Mq 
and 120 Mq as the lower and upper mass limit, respectively. 
The total number of models was 1000 and each model con- 
sisted of 7500 stars which, on average, results in 250 stars 
with m > 10 Mq. We then fitted a power-law MF to each 
of the models in the same way we fitted the observed MF, 
using only the ~ 250 massive stars. In Fig. [S] we show the 
distribution of fitted F- values and from this distribution we 
derive F — — 1.34±0.12. The F-value derived from the obser- 
vations is indicated by the shaded box. Based on this, only 
a small fraction of our models (nine) are consistent with the 
observed MF. However, given the uncertainties in deriving 



4 Harfst et al. 




lower mass limit [M^^^„] 



Figure 5. The distribution of fitted F-values of model MFs. 
Shown are the resuhs from the models (blue pluses) together 
with a fitted Gaussian distribution. The observed MP slope is 
indicated by the shaded area. 



Figure 6. The initial mass of the cluster as a function of the 
lower mass-limit of the IMF for different IMFs. The three lines 
show different normalizations. 



the initial masses of stars (which depend very much on the 
model for the rather unknown mass loss rate) and since so 
far any effects from the dynamical evolution of the cluster 
are not taken into account, we decided to use the slope of 
the IMF as a free parameter and studied, how a different 
IMF affects the present-day MF in the dynamically evolved 
and mass-segregated cluster. We also varied the lower mass 
limit of the IMF to test for a possible truncation of the IMF 
at lower masses. 



3.2 Initial mass of the cluster 

The choice of an IMF is also important for determining 
the total mass of the cluster . Th e present day ma ss of the 
cluster is ^ 1 - 2 ■ 10'^ M p, jFiger et al.. 200 2: St ohe et al.1 
l2008l : lEspinoza et al.ll2009l) . lKim et al.l l|2000[ ) estimate that 
Arches could have lost about half its initial mass due to its 
dynamical evolution in the Galactic center, whic h would give 
an int ial mass of ~ 3T0* Mq. On the other hand. lFiger et all 
l|2002h calculated an upper mass limit of 7 ■ 10* M© for the 
cluster, based on the virial theorem and the observed veloc- 
ity dispersion. 

In Fig. [6] we show the initial cluster mass we derive 
for different IMFs. For each IMF we used three different 
normalizations such that the total number A*'ms of star more 
massive than 20 M0 is 200, 250, and 300. We then varied the 
lower mass lim it of the IMF as i t may be truncated in the 
Arches cluster l|Stohe et al.ll2005^ . 

We find, that the total cluster mass can be initially as 
high as ~ 8 • lO" Mq (Salpeter IMF) or even higher if stars 
below 0.5 M0 have formed in the cluster. For a flat IMF 
(r = —0.9), the total cluster mass does not depend much on 
the lower cut-off and is for Nms = 200 about ^ 2 -J^ Mq 
in good agreement with the findings of iKim et al.l (|2000l l. 



3.3 Size of the cluster 

The Arches is a very dense cluster with a centra l mass den- 
sity of ~ 2 • 10^ Mq pc"^ lEspinoza et al.ll2009l ). The tidal 
radius of the Arches cluster is ~ 1 pc and we determine 



the core radius with r-corc = 0.25 pc using all stars and 
with Tcore = 0.14 pc for stars with m > 10 M©. Following 
ICasertano fc HutI (fl98Sl . the core radius is defined through- 
out this paper as the density-weighted average distance of 
stars to the density center. 

The current concentration of the Arches cluster may be 
the result of its dynamical evolution as the cluster is prob- 
ably evolving towards core collapse. We therefore use the 
initial concentration, parameterized in the King model con- 
centration parameter wo, and size of the cluster, namely its 
virial radius, as two more free parameter. The virial radius 
Tvir is defined as 



1 GM^ 

2 \U\ 



(1) 



with the gravitational constant G, the mass M and poten- 
tial energy U of the cluster. Since U cannot be observed, 
it may be more practical to know that the virial radius is 
proportional to the half mass radius. 



3.4 Summary of model parameters 

To summarize, the initial conditions for our cluster models 
are set by the virial radius, slope and low-mass cutoff of the 
IMF, concentration, and number of massive stars, which are 
varied systematically (see Tab.[T]below and Sec. 4.1). 



4 SETUP AND RESULTS OF SIMULATIONS 

Our aim is to find a numerical model of the Arches clus- 
ter that can explain the observed properties of the cluster. 
In this paper, we neglect the effects of the Galactic poten- 
tial as well as stellar evolution. This allows us to perform 
a large number of simulations in order to constrain some 
of the parameters of our models. The simulations were car- 
ried o ut in the Starlab environm ent using the integrator 
kira l|Portegies Zwart et al.l 1200 ll ). GPUs, graphical pro- 
cessing units, were used to accelerate th e calculations via 
the Sapporo library (|Gaburov et al.|[2009l ). 
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Table 1. List of models. 



Model 


Wo 


IMF 




-^^c luster 


-^cluster 


iv(m > 10 Mq) 




parameter 
















^?vir[pc] 


A''ms 


IKW03F05 


3 


flat 


0.5 


22.9 


6.9 


423 


0.1 


- 1.0 


100 - 300 


IKW03F10 


3 


flat 


1.0 


20.5 


3.7 


421 


0.1 


- 1.0 


100 - 300 


IKW03S05 


3 


Salpeter 


0.5 


52.7 


31.9 


552 


0.1 


- 1.0 


100 - 300 


IKW03S10 


3 


Salpeter 


1.0 


39.7 


12.5 


552 


0.1 


- 1.0 


100 - 300 


IKW03S40 


3 


Salpeter 


4.0 


20.6 


1.9 


540 


0.1 


- 1.0 


100 - 300 


IKW05F05 


5 


flat 


0.5 


22.7 


6.9 


413 


0.1 


- 1.0 


100 - 300 


IKW05F10 


5 


flat 


1.0 


20.2 


3.7 


413 


0.1 


- 1.0 


100 - 300 


IKW05S10 


5 


Salpeter 


1.0 


39.0 


12.5 


545 


0.1 


- 1.0 


100 - 300 


IKW05S40 


5 


Salpeter 


4.0 


20.8 


1.9 


543 


0.1 


- 1.0 


100 - 300 


IKW07F05 


7 


flat 


0.5 


22.9 


6.9 


422 


0.1 


- 1.0 


100 - 300 


IKW07F10 


7 


flat 


1.0 


20.7 


3.7 


421 


0.1 


- 1.0 


100 - 300 


IKW07S10 


7 


Salpeter 


1.0 


39.2 


12.5 


537 


0.1 


- 1.0 


100 - 300 


IKW07S40 


7 


Salpeter 


4.0 


20.8 


1.9 


551 


0.1 


- 1.0 


100 - 300 



Columns are: 1) model name; 2) the dimensionless King concentration parameter Wo 3) the IMF used, where flat and Salpeter refer to 
power-law IMFs with F = —0.9 and F = —1.35, respectively; 4) lower IMF mass limit in Mq; 5) total cluster mass in 10^ Mq; 6) total 
number of stars in the cluster in 10''; 7) total number of massive stars (m > 10 Mq); 8) additional model parameter and range 



4.1 Initial conditions 



We use King models (|Kingj|l96^ in virial equilibrium as our 
initial cluster model. As the initial concentration of the clus- 
ter is not known, it is possible that the cluster has evolved 
to its current compactness from a less concentrated model. 
We therefore decided to use three different values for the 
initial concentration parameter PFo- For each concentration 
parameter we also tested both the Salpeter IMF and an IMF 
with a flat slope of F = —0.9. Furthermore, we also varied 
the lower mass limit of the IMF between 0.5 M© and 4M0. 
In total, we used 13 different sets of initial conditions for the 
cluster (see Tab.[TJ. 

In addition, we also varied two free parameters for each 
of these models: the intial virial radius as rvir and the ini- 
tial number A*'ms of stars with m > 20 M©. The latter is 
used to normalize the total cluster mass and we used five 
different values for A^ms between ICQ and 300 in steps of 50. 
This range covers the 127 stars with m > 20 Mp) found in 
the N ACO data and the ~ 200 stars reported bv lFiger et al.l 
l|2002l l, also taking into account that a significant fraction 
of masssive stars is no longer bound to the cluster or not 
observed. The differences between the two data sets can be 
explained by a different field-of-view (smaller for NACO) 
and the applied selection criteria to determine cluster mem- 
bership. Varying A^ms increases the total number of models 
tested to 65. 

Our TV-body simulations are scale-free A^-body units 
because we neglect stellar evolution and the tidal field. 
In A-body units, the gravitational constant G, the total 
mass of the system, and the virial radius are all set to 
unity l|Heggie fc Mathieul [l98^ . The connection between 
the scale-free A'-body units and physical units is then given 

by 



: 0.0045- 



pc^ 



MoMyr^ 



(2) 



where Ui, Um, and Ut are the A-body units for length, mass, 
and time, respectively. The mass unit Um is naturally given 
by the total mass of the Arches cluster MArchcs- A choice of 



rvir defines Ui, which in turn determines Ut via Eq. [5] The 
age of the cluster is f = 2.5 Myr or according to Eq.[2] 



2.5Myr _ 



0.0045- 



.Ma 



M(T 



— (3) 
pc ' 



in the dimensionless A-body time unit. 

The virial radius was varied between 0.1 and Ipc in 
steps of 0.05 pc. However, instead of increasing the number 
of models by another factor of 19, we make use of the scale- 
free nature of our simulations. For each value of rvir, the 
cluster has to be evolved to a different A^-body time unit 
given by Eq. [3]to reach an age of 2.5 Myr in physical units. 
The A'-body times to match the desired values of rvir can 
be computed at the start of the simulation, and then snap- 
shots are written at these A-body times during this simu- 
lations. Each of the snapshots corresponds to a snapshot of 
the Arches cluster at t = 2.5 Myr with a different rvir- 

Each of the 65 models was created ten times with a 
different random realization, so that a total of 650 simula- 
tions were performed. However, once the simulations were 
finished, we averaged the global properties of each model 
before comparing them with the observations. We also took 
into account the incompleteness of the observations by ran- 
domly removing a few stars from our models a ccording to 
the incompleteness tests by IStolte et al.l l|2005h . However, 
this has only little effect on the overall results as we already 
constrained our data sample to stars that are almost com- 
plete in the observations. In the end, we compared more 
than 1,200 simulation snapshots with the observations. 



4.2 Comparing models and observation 

In order to compare the simulation snapshots with the 
NACO data, we first computed a cumulative mass profile 
for each parameter set. In this process, we averaged the 
ten different random realizations and obtained a single mass 
profile. We only selected stars with m > 10 M© and within 
r < 0.4 pc and we also normalized the profile to the to- 
tal mass inside this radius. Two of the resulting profiles are 
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Figure 7. Comparing the normalized cumulative mass profile of model and observation. The red line shows the model data, with error 
bars indicating the standard deviation from averaging over ten random realization of the same model. The black line the observation. 
Two different models are shown. 
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Figure 8. Comparing the cumulative mass function of model and observation. The red line shows the model data, with error bars 
indicating the standard deviation from averaging over ten random realization of the same model. The green line shows the original IMF 
that was used to setup the ICs for the model. Two different models are shown. 



sliown in Fig. [7] (red line witli error bars) in comparison with 
the observed profile (black line). 

In the following we define a number of fitness param- 
eters / which describe the quality of the fit of the model 
to the observation. These parameters are defined such that 
the values range from zero to unity. A value close to unity 
describes a good fit. 

The mass profile fitness parameter fp is used to quantify 
how well the cumulative mass profile of the mode l fits the ob- 
serva tions. We employ a two-sample KS test fsee lPress et al.l 
1 19921 ) for this comparison. The KS-test returns the probabil- 
ity p, that two data sets are drawn from the same distribu- 
tion. In Fig. [71 the p- value is given for two different models 
and it can be seen that a good fit results in a high p- value 
as expected. In the following, we will use 



(4) 



where a value of fp close to unity describes an excellent fit. 

We also compare the total mass of stars with m > 
10 M0 inside of 0.4 pc. We define the mass fitness param- 



eter 



/m = 1 - 



1 - 



A/mocl 



Ma 



(5) 



as an estimator for how well the observed total mass of the 
cluster is reproduced in the model, /a/ = 1 identifies a per- 
fect match and becomes smaller the more Afmodci deviates 

from MArchcs- 

In Fig. |8l we show the comparison of the observed cu- 
mulative MF for the Arches cluster with our models. The 
integrated MF of the model (red dots) is fitted by a power- 
law MF (dotted blue line) to determine the slope F. Two dif- 
ferent models are shown, which started out with a Salpeter 
IMF (left panel) and a flat IMF (right panel). The IMF used 
for the initial conditions of the model is also plotted in each 
panel (dashed green line) . Because we only take stars within 
0.4 pc into account, the measured (observed) MF can differ 
significantly from the original IMF. Similar to we define 
the MF fitness parameter 



,flMF = 1 — 



1 - 



(6) 
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Figure 9. Quality of fit to observations for varying model parameters rvir and Nms- Data from model IKW03S10 (isolated King model 
with Wo = 3 and a Salpeter IMF with rriiow = 0.5 Mq). Panels show from top left to bottom right the fit to the cluster mass, cumulative 
mass profile, IMF and product of all three fits. A fitness parameter close or equal to unity indicates the best fit models. 



with FArchcs — —0.9. 

The quality of the fit for model IKW03S05 is presented 
in Fig. |9l which shows how fp, Jm, and /imp vary for the 
model parameters rvir and A^ms- The cluster mass is best 
fitted with A^ms « 150 and almost independent of rvir- Only 
models with rvir ^ 0.8 require larger values of A^ms- This 
can be explained as follows: models with small rvir evolve 
faster dynamically and these models have gone through core 
collapse already at an age of 2.5 Myr. As a result, the initial 
dependence of A^ms on rvir , which can still be seen for larger 
values of rvir, is lost. 

The top right panel in Fig. |5] shows which models have 
the best fitting mass profiles. In this case, the best mod- 
els are found within a narrow range of rvir-values with no 
dependence on A^ms- Models with rvir 0.5 pc are, in com- 
parison with the observations, too concentrated (see also 
Fig. [Tl whereas rvir > 0.8 pc results in a too shallow mass 
profile. 

The slope of the IMF depends not as strongly on the 
model parameters (bottom left panel) though the general 
trend is, that the best fits are found for larger values of rvir. 
In the last panel, we show the combined fitness /^n which is 
defined as 

/all = fp ■ fM ■ /iMF. (7) 

From this, we find the best fitting model in the model series 
IKW03S05, which has the virial radius rvir = 0.70 pc and 
A'^MS = 150. The number of stars with m > 10 M© and inside 
R < 0.4 pc in this particular model is 208±12 (averaged over 



ten random realizations), compared to 234 in the observed 
sample. 

In Fig. (TDj we compare /aii for a number of different 
model series. The two panels at the top show the results 
for the models IKW03S10 and IKW03S40. The only diflter- 
ence between these two models and also to model IKW03S05 
shown in Fig. |5]is the lower mass limit of the IMF, which is 
1.0, 4.0, and 0.5 M0, respectively. In each of the three mod- 
els the best fit is in a similar regime of the parameter space, 
with rvir ~ 0.6 — 0.7 pc and A^ms ~ 150. In the other four 
models shown, the best fitting models also lie within a small 
area of the parameter space. 

The two bottom panels show results from models start- 
ing initially with a flat IMF. These models produce not as 
good a fit as models starting with a Salpeter IMF. The rea- 
son is that the slope of the IMF is always flattened by the 
dynamical evolution of the cluster (see also Fig. [S}, so that 
models starting with the observed slope cannot produce a 
good fit. This result is independent of the initial concentra- 
tion of the cluster though more concentrated clusters give 
slightly better results. 

We have summarized the results in Tab. [2] where for 
each model series the best values for rvir and A^ms are given 
together with the corresponding value of /aii. We also re- 
peated the whole analysis using all stars down to 4M0 in 
the comparison. In both cases, models with the highest val- 
ues of /all all have a Salpeter IMF. These best fit models 
also have similar value of rvir and A^ms- We computed the 
average of these values weighted with /aii which results in 
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Figure 10. Quality of fit to observations for different models and varying model parameters rvir and -/Vms- The combined fitness f^n is 
shown for models IKW03S10, fKW03S40, fKW03F05, and fKW07F10 from top left to bottom right, respectively 



Table 2. Best fitting models. 



Model 


parameter 


/all 






Rvii [pc] 




m > 10 M© m 


> 4M0 


IKW03F05 


0.65 


200 


0.37 


0.45 


IKW03F10 


0.7 


200 


0.43 


0.47 


IKW03S05 


0.7 


150 


0.85 


0.68 


IKW03S10 


0.65 


150 


0.79 


0.82 


IKW03S40 


0.6 


200 


0.66 


0.56 


IKW05F05 


0.85 


200 


0.46 


0.45 


IKW05F10 


0.9 


250 


0.48 


0.42 


IKW05S10 


0.8 


200 


0.69 


0.64 


IKW05S40 


0.75 


200 


0.88 


0.77 


IKW07F05 


0.75 


200 


0.52 


0.43 


IKW07F10 


0.85 


250 


0.62 


0.41 


IKW07S10 


1.0 


200 


0.66 


0.76 


IKW07S40 


0.85 


250 


0.78 


0.70 



Columns are: 1) model name; 2) values of additional model pa- 
rameter; 3) and 4) /^n comparing stars with m > 10 Mq and 

m > 4M0 



rvir = 0.76 ± 0.12 pc and A^ms = 192 ± 32 for the compar- 
ison using stars with m > 10 M© and rvir = 0.62 ± 0.15 pc 
and A^MS = 164 ± 35 for stars with m > 4Mq. The ini- 
tial mass of the cluster depends on the lower mass limit of 
the IMF which is not constrained by our models unfortu- 
nately. Based on our results for stars with m > 10 Mq, we 
get M = (51 ±8) • lO^^M©, M = (38 ± 6) ■ 10^ Mq, and 



M = (20 ± 3) • 10^ Mq for lower mass limits of the IMF of 
0.5, 1, 4M0. 



5 DISCUSSION 

A number of models in Tab. [2] can be considered best fit 
models, for example model IKW03S05 with rvir = 0.70 and 
A'^MS ~ 150. A snapshot of this model at t = 2.5 Myr is 
shown in Fig. 1111 Stellar masses are indicated by the point 
size and gray-scale and the center of density of the cluster is 
located at the origin. The dashed circle has a radius of 0.4 pc. 
In the left panel, the NACO data is plotted for comparision 
in the same way as the simulation data in the right panel. 
In the latter, stars have been randomly removed to mimic 
incompleteness. The probability that a star is removed is 
given by a function fitted to the data shown in Fig. [2] and 
depends on mass and position of the star. 

At 2.5 Myr, the cluster has a total mass of - 2 ■ 10" M0 
inside a radius of 0.4 pc and twice that mass inside the tidal 
radius of 1 pc. On average, our best models predict that 
the Arches cluster is more massive than observed, even if 
we take into account mass loss by stellar evolution. About 

6 ■ 10^ Mq of the total mass are in stars below 1 Mq . Our 
results are in good ag r eemen t with the previous findings of 
iPortegies Zwart et al.i (|2002l '). 

Generally, only models with the Salpeter IMF produce 
acceptable fits and we therefore conclude that the IMF in the 
Arches cluster is normal despite the extreme environment in 
which the cluster is forme d. However, we can not rule out a 
turn-over as suggested by jStolte et al.ll200"5h as our results 
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Figure 11. The observed cluster (left) and a snapshot of one of the best fitting models (right) in comparison. The images are centered 
on the center of density and the circles indicate a radius of 0.4 pc. Gray-scale and point size represent stellar masses. In the right panel, 
stars have been removed randomly to mimic incompleteness. 
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Figure 12. The IMF of the best fitting model for all stars and 
the stars in two radial bins. 



are not very sensitive to the lower mass limit of the IMF. 
Our best fitting models can have any lower mass limit in the 
the range of 0.5 — 4.0 Mq that we investigated. 

Fig. ll2l shows the MF of one of our best fit models. A bin 
size of Alog{m/M.Q) = 0.2 was used to obtain the individual 
star counts and the data was also averaged over the ten dif- 
ferent realization of this model. The top line shows the MF 
for all stars which is simply the Salpeter IMF used in the 
initial conditions of this model. The other tw o lines show 
the M F in two radial bins (the s ame used by IStolte et al.l 
l|2005l ) and lEspinoza et all (120091 ')'). The MF is noticeably 
flattened for massive stars (m > 10 M©) in the inner radial 
bin {R < 0.2 pc), while it remains unchanged in the outer 
radial bin (0.2 < R < 0.4 pc). The flattening of the MF in 
the cluster core has been reported in all the previous obser- 
vational studies, though with varying F-values determined 
for the slope. Our best flt models are in v ery good agree- 
ment, however, with the latest results from lEspinoza et al.l 
l|20090 . Stars with m > 10 M© can be found in our models as 
far as 10 pc from the cluster center. However, we find that in 
order to measure the correct slope of the MF, we only need 
the stars inside the tidal radius of 1 pc. 
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Figure 13. The core radius evolution with time of model 
IKW03S05 with r^ir = 0.7pc and A^ms = 150. The black line 
with errorbars shows the average of ten random realizations of 
the same model (grey lines) . Core collapse happens at ~ 4 Myr. 



No turn-over at the low mass end of the IMF can be 
seen. This means that the turn-over seen by IStolte et al.l 
(j2005l ) cannot be explained by mass segregation (unlike the 
flat IMF at the high-mass end). Deriving the IMF from our 
simulation data is not hampered by incompleteness, selec- 
tion effects, and mass determination from observed lumi- 
nosities. All these effects may bias the observed star counts 
to produce a turn-over, however, the turn-over appears at 
~ 6 M0 where the data is still 50% complete. So, either the 
IMF is indeed truncated or another effect has to be consid- 
ered. One possible explanation is that tidal stripping prefer- 
ably removes low-mass stars from the cluster. This effect is 
not included in our current simulati ons but will be tested 
in a follow-up paper. Alternatively, lEspinoza^ et all (|2009l ) 
pointed out that local variations in extinction can account 
for (some of) the flattening of the obserevd MF. The same 
effect could cause lower-mass stars to be lost preferentially. 
Then the completeness function, which only takes into ac- 
count crowding and sensitivity effect, would not be sufficient 
to correct the low-mass end of the IMF. 
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In Fig. [13] we show the core radius evolution with time 
for one of our best fit models. At its current age of 3.5 Myr 
the cluster is a little more than halfway to core collapse 
which will occur at ~ 4 Myr. This result is in agreement 
with the findings of IPortegies Zwart et al.l l|2007h . 



6 SUMMARY 

We have performed a large number of A'^-body simulations 
in order to find the best fitting model for the Arches cluster. 
The available observational data has been used to constrain 
the free parameters in our model. In a systematic analysis, 
we compared the total mass, the cumulative mass profile, 
and the present-day MF and defined fitness parameters for 
each of the three observables. 

The main conclusion from our analysis is that the 
Arches cluster, despite of being born in an extreme envi- 
ronment, has formed with a standard Salpeter IMF. Due to 
mass segregation, the slope of the observed MF is flattened 
inside a radius of 0.4 pc. This radius was imposed a limit 
from the obversavational data used, and we estimate that a 
limiting radius of ~ 1 pc would be required for the observed 
MF to match the underlying IMF. 

We conclude from our best fitting models that the 
Arches cluster was born with an initial virial radius of 0.76± 
0.12 pc and an initial total mass of about (5.1±0.8)-10'' M© 
assuming a lower mass limit of 0.5 M© for the Salpeter IMF. 
The lower mass limit cannot be constrained well from our 
models, giving rise to additional uncertainties in determin- 
ing the initial cluster mass. The King model concentration 
parameter of the best fitting models is Wo = 3 — 5, how- 
ever, the Wo = 7 produced also reasonable fits so that this 
parameter is not very well constrained. 

We neglected the Galactic potential and also stellar evo- 
lution for the simulations in this paper. These processes will 
be included in a following paper to get a more realistic model 
of the Arches cluster, however this first step was needed 
in order to reduce the number of free parameters for the 
these (computationally more expensive) simulations, which 
further constrain the dynamical evolution and tidal effects 
acting on the Arches cluster, and thereby the initial condi- 
tions of this nearby starburst, such as the cluster mass, the 
orbital motion and the IMF at the birth of the cluster. 
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